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We present a group theoretical study of the symmetry- 
broken unrestricted Hartree-Fock orbitals and electron den- 
sities in the case of a two-dimensional A''-electron single 
quantum dot (with and without an external magnetic field). 
The breaking of rotational symmetry results in canonical 
orbitals that (1) are associated with the eigenvectors of 
a Hiickel hamiltonian having sites at the positions deter- 
mined by the equilibrium molecular configuration of the clas- 
sical A''-electron problem, and (2) transform according to 
the irreducible representations of the point group specified 
by the discrete symmetries of this classical molecular con- 
figuration. Through restoration of the total-spin and ro- 
tational symmetries via post-Hartree-Fock projection tech- 
niques, we show that the point-group discrete symmetry 
of the unrestricted Hartree-Fock wave function underlies 
the appearance of magic angular momenta (familiar from 
exact-diagonalization studies) in the excitation spectra of 
the quantum dot. Furthermore, this two-step symmetry- 
breaking/symmetry-restoration method accurately describes 
the energy spectra associated with the magic angular mo- 
menta. 



Pacs Numbers: 73, 73.21-b, 73.21.La 



I. INTRODUCTION 

A. Background on the mean-field breaking of spatial 
symmetries in quantum dots 

Two-dimensional (2D) quantum dots (QD's) created 
at semiconductor interfaces with refined control of their 
size, shape and number of electrons are often referred^"^ 
to as "artificial atoms". For high magnetic fields {B), 
it has been known for some time (ever since the pio- 
neering work^ of Laughlin in 1983 concerning the frac- 
tional quantum Hall effect) that 2D few-electron systems 
exhibit complex strongly-correlated many-body physics. 
Nevertheless, for low magnetic fields, the term artificial 
atoms was used initially to suggest that the physics of 
electrons in such man-made nanostructures is exclusively 
related to that underlying the traditional description^ of 
natural atoms (pertaining particularly to electronic shells 
and the Aufbau principle), where the electrons are taken^ 
to be moving in a spherically averaged effective central 
mean field. This traditional picture was given support 
by experimental studies'^'^ on vertical QD's, which were 



followed^ by a series of sophisticated theoretical investi- 
gations yielding results conforming to it. 

However, in 1999, the circular (for 2D QD's) central- 
mean-field picture was challenged by the discovery of 
solutions with broken space symmetry in the context 
of spin-and-space unrestricted Hartree-Fock (sS-UHF) 
mean-field calculations. These broken symmetry (BS) 
solutions appear spontaneously (due to a phase transi- 
tion) when the strength of the interelectron repulsion 
relative to the zero-point kinetic energy (Rw) exceeds a 
certain critical value. They indicate formation of Wigner 
(or electron) molecules (WM's or EM's) with the elec- 
trons located at the vertices of nested regular polygons 
(often referred to as concentric rings), familiar from stud- 
ies of classical point charges. ^^'^^ Such molecules were 
characterized by us as strong or weak Wigner crystal- 
lites depending on their rigidity or lack thereof (i.e., flop- 
piness), suggesting the existence of additional^'^ "phase 
transitions" as a function of the parameter Rw (see sec- 
tion H.B for its precise definition). Furthermore, it was 
noted^ that the symmetry breaking should be accompa- 
nied by the emergence of a spectrum of collective rovi- 
brational excitations (finite-size analogs of the Goldstone 
modes) . A subsequent investigation^* based on exact so- 
lutions for a Helium QD (2e QD) confirmed these results 
and provided a systematic study of the molecular rovi- 
brational collective spectra and their transition to inde- 
pendent particle excitations, as the rigidity of the WM 
was reduced through variation of the controlling param- 
eter Rw- 

We remark here that the lower-energy BS UHF so- 
lutions already capture^° part of the correlation en- 
ergy, compared to the restricted HF (RHF) ones. Im- 
proved numerical accuracy has been achieved in subse- 
quent studies^^""'^^ through the restoration of the broken 
symmetries via projection techniques. Consequently, the 
methodology of symmetry breaking at the UHF mean- 
field level and of subsequent symmetry restoration via 
post Hartree-Fock methods^* provides a systematic con- 
trolled hierarchy of approximations toward the exact so- 
lution, with anticipated advantages for complex many 
electron systems (under field-free conditions'^ and in the 
presence of a magnetic field'^'^^), whose treatment is 
computationally prohibitive with other methods (e.g., ex- 
act diagonalization). 
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B. Background on the mean- field breaking of 
symmetries in other finite-size fermionic systems 

The mean field approach has been a useful tool in eluci- 
dating the physics of small fermionic systems, from natu- 
ral atoms and atomic nuclei to metallic nanoclusters, and 
most recently of two-dimensional quantum dots. Of par- 
ticular interest for motivating the present work (due to 
spatial-symmetry-breaking aspects) has been the mean- 
field description of deformed nuclei and metal clusters 
(exhibiting ellipsoidal shapes). At a first level, defor- 
mation effects in these systems can be investigated via 
semi-empirical mean-field models, like the particle-rotor 
model^^ of Bohr and Mottelson (nuclei), the anisotropic- 
harmonic-oscillator model of Nilsson (nuclei^" and metal 
clusters^^), and the shell-correction method of Strutin- 
sky (nuclei^^ and metal clusters^^'^^). At the microscopic 
level, the mean field is often described^'"''^^ via the self- 
consistent single-determinantal Hartree-Fock (HF) the- 
ory. At this level, however, the description of defor- 
mation effects mentioned above requires^^'^'' consider- 
ation of unrestricted Hartree-Fock wave functions that 
break explicitly the rotational symmetries of the origi- 
nal many-body hamiltonian, but yield HF determinants 
with lower energy compared to the symmetry-adapted 
restricted Hartree-Fock solutions. 

In earlier publications, we have shown that, in 
the strongly correlated regime, UHF solutions violating 
the rotational (circular) symmetry arise most naturally 
in the case of 2D single QD's and for both the cases of 
zero and high magnetic field. Unlike the case of atomic 
nuclei, however, where symmetry breaking is associated 
with quadrupole deformations, spontaneous symmetry 
breaking in 2D QD's induces electron localization asso- 
ciated with formation of WM's. 

The violation in the mean-field approximation of the 
symmetries of the original many-body hamiltonian ap- 
pears to be paradoxical at a first glance. However, for the 
specific cases arising in Nuclear Physics and Quantum 
Chemistry, two theoretical developments have resolved 
this paradox. They are: (1) the theory of restoration 
of broken symmetries via projection techniques, ^^'"^^ and 
(2) the group theoretical analysis of symmetry-broken HF 
orbitals and solutions in chemical reactions initiated by 
Fukutome and coworkers, who used of course the sym- 
metry groups associated with the natural 3D molecules. 
Despite the different fields, the general principles estab- 
lished in these earlier theoretical developments have pro- 
vided a wellspring of assistance in our investigations of 
symmetry breaking in QD's. In particular, the restora- 
tion of broken circular symmetry in the case of single 
QD's has already been demonstrated by us in three re- 
cent publications. ^^"^'''^^ 



C. Content of this paper 

In the present paper, we will provide an in-depth 

group theoretical analysis of broken-symmetry UHF or- 
bitals and electron densities (ED's) in the case of sin- 
gle parabolic QD's. We will show that such an anal- 
ysis provides further support for our earlier interpreta- 
tion concerning the spontaneous formation of collectively 
rotating electron (or Wigner) molecules (REM's). In- 
deed we will demonstrate deep analogies between the 
electronic structure of the WM and that of the nat- 
ural 3D molecules. In particular, we will show that 
the breaking of rotational symmetry results in canon- 
ical UHF orbitals that are associated with the eigen- 
vectors of a molecular-type Hiickel hamiltonian having 
"atomic" sites at positions specified by the equilibrium 
configuration of the classical A''-electron problem; these 
"atomic" sites are localization sites for the electrons, and 
they do not imply the presence of a positive nucleus. 
Thus, in contrast to the fully delocalized and symmetry- 
adapted, independent-particle-model-type orbitals of the 
RHF, the BS UHF orbitals with the same spin direction 
resemble closely the molecular orbitals (MO's) formed by 
linear combinations of atomic orbitals (LCAO's), which 
are prevalent^^ in Chemistry. (Naturally, the LCAO be- 
havior of the UHF orbitals with the same spin direc- 
tion allows for a more precise understanding of the term 
"electron localization" used by us in previous publica- 
tions.) An important conclusion of the present paper 
is that the BS UHF orbitals are not necessarily unique; 
what matters, in analogy with the LCAO-MO's of nat- 
ural molecules, is that they transform according to the 
irreducible representations of the point group specified 
by the discrete symmetries of the classical equilibrium 
configuration.^^ 

Our study leads to the following two main results: (i) 
in analogy with 3D natural molecules, the WM's can ro- 
tate and the restoration of the total-spin and rotational 
symmetries via projection techniques describes their low- 
est rotational bands {yrast bands^'*'"^"') in addition to 
the ground state, and (ii) the lowering of the symme- 
try, which results in the (discrete) point-group symme- 
try of the UHF wave function, underlies the appearance 
of the sequences of magic angular momenta (familiar 
from exact-diagonalization studies'^"'' *'^) in the excitation 
spectra of single QD's. Since exact-diagonalization meth- 
ods are typically restricted to small sizes with N < 10, 
the present two-step method of breakage and subsequent 
restoration of symmetries offers a promising new av- 
enue for accurately describing larger 2D electronic sys- 
tems. A concrete example of the potential of this ap- 
proach is provided by Ref. 16 and Ref. 17, where our 
use of the the symmetry-breaking/symmetry-restoration 
method yielded analytic expressions for correlated wave 
functions that offer a better description of the A/'-electron 
problem in high magnetic fields compared to the Jastrow- 
Laughlin^ expression. 



2 



Since the group-theoretical aspects of symmetry break- 
ing at the mean-field level (and their relation to the prop- 
erties of the exact solutions) remain a vastly unexplored 
territory in the area of condensed-matter finite-size sys- 
tems, in the following we will present an introductory 
investigation of them through a series of rather simple, 
but nontrivial, illustrative examples from the field of 2D 
parabolic QD's. The plan of the paper is as follows: Sec- 
tion II reviews briefly the set of UHF equations employed 
by us; Section III and section IV present the case of three 
electrons in the absence and in the presence of an exter- 
nal magnetic field B, respectively. The more complicated 
case of six electrons at B = is investigated in section 
V, while section VI discusses the companion step of the 
restoration of broken symmetries, which underlies the ap- 
pearance of magic angular momenta in the exact spectra. 
Finally, section VII presents a summary. 

Before leaving the Introduction, we wish to stress that 
analogs of several of the group theoretical concepts and 
manipulations employed in this paper can be found in 
textbooks concerning standard applications of symmetry 
groups to the electronic structure and chemistry of 3D 
natural molecules. Here, however, we will use these oth- 
erwise well-known group theoretical aspects to elucidate 
the molecular interpretation of the BS UHF determinants 
(and associated orbitals) in the unexpected context of 
a newly arising area of physics, namely the physics of 
strong correlations in 2D circular artificial atoms (QD's). 

II. THE UHF EQUATIONS 

A. The Pople-Nesbet equations 

The UHF equations we are using are the Pople- 
Nesbet^^ equations described in detail in Ch. 3.8 of Ref. 

26. For completeness, we present here a brief descrip- 
tion of them, along with details of their computational 
implementation by us to the 2D case of semiconductor 
QD's. 

The key point is that electrons of a (up) spin are de- 
scribed by one set of spatial orbitals = 1) 2, K}, 
while electrons of (3 (down) spin are described by a differ- 
ent set of spatial orbitals {ipj\j = 1,2,...,K} (of course 
in the RHF t/;" — i^j = i'j)- Next, one introduces a set 
of basis functions {VjuIm = 1,2,..., if} (constructed to 
be orthonormal in our 2D case), and expands the UHF 
orbitals as 

K 

^T = J2'^>^^ i = l,2,...,i^ (1) 

K 

V'f = EO'^' i = h2,...,K. (2) 



The UHF equations are a system of two coupled matrix 
eigenvalue problems, 

F"C" = C"E" (3) 

F^C<^ = C^E^, (4) 

where F"^'^) are the Fock-operator matrices and C"*-^-* 
are the vectors formed with the coefficients in the ex- 
pansions (1) and (2). The matrices E"^''^ are diagonal, 
and as a result equations (3) and (4) are canonical (stan- 
dard). Notice that noncanonical forms of HF equations 
are also possible (see Ch. 3.2.2 of Ref. 26). Since the 
selfconsistent iterative solution of the HF equations can 
be computationally implemented only in their canonical 
form, heretofore canonical orbitals and solutions will al- 
ways be implied, unless otherwise noted explicitly. We 
note that the coupling between the two UHF equations 
(3) and (4) is given explicitly in the expressions for the 
elements of the Fock matrices below [Eqs. (7) and (8)] . 

Introducing the density matrices P"(^) for q:(/3) elec- 
trons, 

p'^. = Y.^';iai.czr (5) 

a 

P^v = X] ^^ai^va)*^ (6) 
a 

where iV" -I- = N, the elements of the Fock-operator 
matrices are given by 

f;, = h^^ + J2Y1 PxAif^^W^) - (m^Ia^^)] 

A cr 
A a- 

A CT 

+ EE^a.(m^I^A), (8) 

A o- 

where H^^ are the elements of the single electron hamil- 
tonian with an external magnetic field B and an appro- 
priate potential confinement, and the Coulomb repulsion 
is expressed via the two-electron integrals 




dridr2V3;(ri)(/5;(r2)T— — — rV'i^(ri)(^A(r2), (9) 



with K being the dielectric constant of the semiconductor 
material. Of course, the Greek indices n, v, A, and a run 
from 1 to i^. 
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FIG. 1. UHF electron density in a parabolic QD for A'' = 19 
and Sz = 19/2, exhibiting breaking of the circular symmetry 
at Rw = 5 and B = 0. The choice of the remaining param- 
eters is: Hljo = 5 meV and m* — 0.067me. Distances are in 
nanometers and the electron density in 10"* nm"'^. 



For a QD, the external confinement is assumed to be 
parabolic, and the single-particle hamiltonian in a per- 
pendicular external magnetic field B is written as 



H = 



(P 



-eA/c)^ 1 * ? 2x 9*f^BT^ 



(10) 



The vector potential A is given in the symmetric gauge 
by 

A{r) = ^Bxr=^{-By,Bx,0), (11) 

and the last term in Eq. (10) is the Zeeman interaction 
with g* being the effective g factor, /z_b the Bohr magne- 
ton, and s the spin of an individual electron, m* is the 
effective electron mass and loq is the frequency parameter 
of the parabolic potential confinement. 

The system of the two coupled UHF matrix equations 
(3) and (4) is solved selfconsistently through iteration 
cycles."*^ For obtaining the numerical solutions, we have 
used a set of i^T = 78 basis states (fi's that are chosen 
to be the product wave functions formed out from the 
eigenstates of the one-dimensional harmonic oscillators 
along the x and y axes. Note that the value K = 78 cor- 
responds to all the states of the associated 2D harmonic 
oscillator up to and including the 12th major shell. 

Having obtained the selfconsistent solution, the total 
UHF energy is calculated as 



(12) 



B. Solutions representing Wigner molecules 

As a typical example of a solution that can be ex- 
tracted from the above UHF equations, we mention the 
case of = 19 electrons for huJo = 5 meV, Rw — 5 
(k = 3.8191), and B = 0. The Wigner parameter 
Rw = Q/S-^Oj where Q is the Coulomb interaction 
strength; Q = e^/ulo, with = {Ti/m*LjQ)^/'^ being the 
spatial extent of the lowest single-electron wave function 
in the parabolic confinement. 

Fig. 1 displays the total electron density of the BS UHF 
solution for these parameters, which exhibits breaking of 
the rotational symmetry. In accordance with ED's for 
smaller dot sizes published by us earlier,^' the ED in 
Fig. 1 is highly suggestive of the formation of a Wigner 
molecule, with an (1,6,12) ring structure in the present 
case. This polygonal ring structure agrees with the clas- 
sical one^^ and is sufficiently complex to instill confidence 
that the Wigner-molecule interpretation is valid. The fol- 
lowing question, however, arises naturally at this point: 
is such molecular interpretation limited to the intuition 
provided by the landscapes of the total ED's, or are there 
deeper analogies with the electronic structure of natural 
3D molecules? The answer to the second part of this 
question is in the positive, and the remainder of this 
paper is devoted to discovering such analogies. How- 
ever, since the TV = 19 case represents a rather com- 
plicated group-theoretical structure, for simplicity and 
transparency, we will study in the following smaller QD 
sizes. This, however, will not result in any loss of gener- 
ality in our conclusions. 

In previous publications, we found that space sym- 
metry breaking in the UHF equations appears sponta- 
neously for Rw > 1- Below we choose to work with larger 
values of Rw (e.g., 10 or higher), for which the effects 
of strong correlations are fully developed. The group- 
theoretical investigation of the intermediate regime near 
the phase transition is left for a future publication. In all 
calculated cases, we used Hujo — 5 meV, m* = 0.067me 
(GaAs), and g* = -0.44 (GaAs). 

We note that the Pople-Nesbet UHF equations are pri- 
marily employed in Quantum Chemistry for studying the 
ground states of open-shell molecules and atoms. Unlike 
our studies of QD's, however, such chemical UHF studies 
consider mainly the breaking of the total spin symmetry, 
and not that of the space symmetries. As a result, for 
purposes of emphasis and clarity, we have often used (see, 
e.g., our previous papers) prefixes to indicate the specific 
unrestrictions involved in our UHF solutions, i.e., the 
prefix s- for the total-spin and the prefix S- for the space 
unrestriction. 
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FIG. 2. UHF solution in a parabolic QD exhibiting a pure 
spin density wave for A'^ = 14, Sz = 0, Rw ~ 0.8, and 
5 = 0. (a) The total electron density exhibiting circular 
symmetry; (b) The spin density exhibiting azimuthal modu- 
lation (note the 12 humps whose number is smaller than the 
number of electrons). The choice of the remaining parame- 
ters is: hiOo = 5 meV and m* = 0.067me. Distances are in 
nanometers and the electron and spin densities in 10~* nm"'^. 

C. Solutions representing pure spin density waves 

Before leaving this section, we mention another class 
of BS solutions, which can appear in single QD's, namely 
the spin density waves (SDW's). The SDW's are unre- 
lated to electron localization and thus are quite distinct 
from the WM's; in single QD's, they were obtained'*^ 
earlier within the framework of spin density functional 
theory."*^ To emphasize the different nature of SDW's 
and WM's, we present in Fig. 2 an example of a SDW 




FIG. 3. The RHF case for = 3 and Sz = 3/2 at 
Rw = 10. (a-c): real orbitals at B = 0. (e-g): the mod- 
ulus square of the complex orbitals for a vanishingly small 
value B — 0.0001 T. (d): the corresponding circular electron 
density for both cases. The choice of the remaining param- 
eters is: huJo = 5 meV and m* = 0.067me. Distances are 
in nanometers. The real orbitals are in 10"'' nm~^ and the 
orbital densities and total ED in 10~* nm~^. The arrows 
indicate the spin direction. 

obtained with the UHF approach [the corresponding pa- 
rameters are: N = 14, 5"^ = 0, = 0.8 (n = 23.8693), 
and B = 0]. Unlike the case of WM's, the SDW exhibits 
a circular ED [see Fig. 2(a)], and thus does not break 
the rotational symmetry. Naturally, in keeping with its 
name, the SDW breaks the total spin symmetry and ex- 
hibits azimuthal modulations in the spin density^^ (SD) 
[see Fig. 2(b); however, the number of humps is smaller 
than the number of electrons^^] . The SDW's in single 
QD's appear for Rw < 1 and are of lesser importance;'''' 
thus in the following we will exclusively study the case of 
WM's. 
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FIG. 4. The S-UHF case exhibiting breaking of the circular 
symmetry for TV = 3 and 5*2 — 3/2 at Rw = 10 and B = 0. 
(a-c): real orbitals. (d): the corresponding electron density. 
The choice of the remaining parameters is: huo = 5 meV 
and m* = 0.067me. Distances are in nanometers. The real 
orbitals are in 10~^ nm~^ and the total ED in 10~* nm~^. 
The arrows indicate the spin direction. 



III. THREE ELECTRONS WITHOUT 
MAGNETIC FIELD 

A. The Sz ~ 3/2 fully spin polarized case 

We begin with the case of = 3 fuUy spin polarized 
{Sz =3/2) electrons in the absence of a magnetic field 
and for Rw = 10 (k = 1.9095). Fully spin polarized UHF 
determinants preserve the total spin, but for this value of 
Rw the lowest in energy UHF solution is one with broken 
circular symmetry. As it will be seen below, broken rota- 
tional symmetry does not imply no space symmetry, but 
a lower point-group symmetry. Before proceeding with 
the study of the BS solution, however, it will be helpful to 
review the symmetry-adapted RHF solution first. This 
RHF solution can be obtained from the UHF equations 
(3) and (4) by using a circular electron-density guess as 
the input of the first iteration. In the independent par- 
ticle model, the N = 3, = 3/2, and S = 2D case 
corresponds to a closed electronic shell with configuration 



Islp+lp- or lslpj:lpy {p± (X Px±iPy 



) , and thus 



the independent-particle-model ED is necessarily circu- 
lar. We will confirm that the RHF solution conforms in- 
deed to the prediction of the independent particle model, 
and subsequently we will contrast the UHF solution to 
it. 

The RHF solution for Rw = 10 and B = has an en- 
ergy of 92.217 meV; the corresponding orbitals are real 
and are displayed in the left column of Fig. 3. They are 



like the Is [Fig. 3(a)], Ip, [Fig. 3(b)], and \py [Fig. 3(c)] 
orbitals of the independent particle model. The node- 
less Is orbital has a maximum at r > due to the large 
Coulomb repulsion; its energy is 44.526 meV. The energy 
of the two degenerate Px and py orbitals is 50.489 meV. 
Notice that neither the px nor py orbital is rotationally 
symmetric; however, the total ED [Fig. 3(d)] has the ex- 
pected circular symmetry. It is of interest to obtain the 
RHF solution for a very small external magnetic field 
(i.e., in the limit i? — > 0). In this case, the calculated 
total and orbital energies, as well as the total ED [Fig. 
3(d)] remain unchanged. However, the two degenerate p 
orbitals are now complex [p±, see Fig. 3(f) and Fig. 3(g)] 
and have good angular momenta Z = ±1, and thus their 
modulus square is circularly symmetric. 

We focus now on the sS-UHF solution for iV = 3 and 
Sz = 3/2, and for the same choice of parameters as with 
the RHF case. The UHF total energy is 89.691 meV, 
and thus it is lower than the corresponding RHF one. In 
Fig. 4 we display the UHF symmetry-violating orbitals 
(a— c) whose energies are (a) 44.801 meV, (b) and (c) 
46.546 meV, namely the two orbitals (b) and (c) with 
the higher energies are again degenerate in energy. An 
inspection of Fig. 4 immediately reveals that these or- 
bitals have retained some properties of the delocalized 
Is, \px and Ipy orbitals, familiar from the independent 
particle model and the RHF; that is, orbital (a) is node- 
less, while each one of the orbitals (b) and (c) has a single 
nodal line. However, overall the BS orbitals (a— c) drasti- 
cally differ from the orbitals of the independent particle 
model. In particular, they are associated with specific 
sites (within the QD) forming an equilateral triangle, and 
thus they can be described as having the structure of a 
linear combination of "atomic" (site) orbitals (LCAO's). 
Such LCAO molecular orbitals (MO's) are familiar in 
natural molecules, and this analogy supports the term 
"electron (or Wigner) molecules" for characterizing the 
BS UHF solutions. 

In the LCAO-MO approximation, one needs to solve 
a matrix eigenvalue equation determined by the overlaps 
^ij j) ^'^d the hamiltonian matrix elements Hij 

and Hii between the atomic orbitals. A further simpli- 
fied approximation,^® the Hiickel approximation (HiiA), 
consists in taking all 5*^ = 0, and all Hij = unless 
the zth and jth atoms (sites) are adjacent. For our ex- 
ample this latter approximation is applicable, since the 
value Rw = 10 is rather high. When using the notations 
e^Hii = H22 = H33 and -/? = Hu = H13 = H23 < 0, 
the Hiickel eigenvalue equation for the case of = 3 elec- 
trons on the vertices of an equilateral triangle is written 
as 



(13) 



and the associated LCAO-MO orbitals are tpi = fl4>i + 
/2^2 + /3^3, having Ei eigenvalues with i — 1, 2, 3. The 






h 


h 


= E\ f2 


h 
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(pj^s arc the original Gaussian-type atomic (site) orbitals. 

From the eigenvalues and eigenvectors of (13), one 
finds the following three LCAO-MO's: 

V-i = (01 + </'2 + <^3)/V3 (14) 

with energy Ei = e — 2/3, 

i>2 = (201 - 02 - 03)/V6 (15) 

with energy £'2 = e + /3, and 

V'3 = (02 - 03)/\/2 (16) 

with energy £"3 = E2. It is apparent that the struc- 
ture of these three LCAO-MO's and the level diagram 
of their energies agrees very well with the corresponding 
symmetry-broken UHF orbitals [displayed in Figs. 4(a) 
— 4(c)] and their energies. (Using the HF values for Ei 
and E2, one finds e « 45.961 meV and /3 « 0.585 meV.) 
We notice here that such LCAO orbitals are familiar in 
Organic Chemistry and arc associated with the theoreti- 
cal description of Carbocyclic Systems, and in particular 
the molecule C3H3 (cyclopropenyl, see, e.g., Ref. 33). 

Naturally, since the orbitals (b) and (c) arc degener- 
ate in energy, they are not uniquely defined: any linear 
combination associated with a unitary 2x2 transforma- 
tion will produce a pair of different, but equivalent (b') 
and (c') orbitals. The fact that the UHF orbitals in Fig. 
4 have the specific highly symmetrized (sec below) form 
given above is the result of an accidental choice of the ini- 
tial electron-density input in the HF iteration. We have 
checked that any such pair of (b') and (c') orbitals leaves 
the 2D total UHF electron density unchanged. This sug- 
gests that there is an underlying group theoretical struc- 
ture that governs the BS UHF orbitals. The important 
point is not the uniqueness or not of the 2D sS-UHF 
orbitals, but the fact that they transform according to 
the irreducible representations of specific point groups, 
leaving both the sS-UHF determinant and the associated 
electron densities invariant. Given the importance of this 
observation, we proceed in the rest of this section with a 
group theoretical analysis of the BS sS-UHF orbitals for 
the TV = 3 and Sz — 3/2 case. 

The ED portrayed in Fig. 4(d) remains invariant under 
certain geometrical symmetry operations, namely those 
of an unmarked, plane and equilateral triangle. They are: 
(1) The identity E; (2) The two rotations C3 (rotation 
by 27r/3) and C| (rotation by 47r/3); and (III) The three 
reflections a^, (tI/ , and a^^^ through the three vertical 
planes, one passing through each vertex of the triangle. 
These symmetry operations for the unmarked equilateral 
triangle constitute the elements of the group Csy.^^'"^^ 

One of the main applications of group theory in Chem- 
istry is the determination of the eigenfunctions of the 
Schrodinger equation without actually solving the matrix 
equation (13). This is achieved by constructing the so- 
called symmetry- adapted linear combinations (SALC's) 



TABLE I. Character table for the cyclic group C3 
[e = exp(27ri/3)] 





E 


Cs 


ci 


A 


1 


1 


1 


E' 


1 


e 


e* 


E" 


1 


e* 


e 



of AO's. A widely used tool for constructing SALC's is 
the projection operator 

where R stands for any one of the symmetry operations 
of the molecule, and x^(i?) are the characters of the /ith 
irreducible representation of the set of R's. (The x'^'s 
are tabulated in the socalled character tables. ^^'^^) \G\ 
denotes the order of the group and the dimension of 
the representation. 

The task of finding the SALC's for a set of three l.s- 
type AO's exhibiting the C^y symmetry of an equilateral 
triangle can be simplified, since the pure rotational sym- 
metry by itself (the rotations C3 and C|, and not the 
refiections a„'s through the vertical planes) is sufficient 
for their determination. Thus one needs to consider the 
simpler character table^° of the cyclic group C3 (see Ta- 
ble I) 

From Tabic I. one sees that the set of the three Is AO's 
situated at the vertices of an equilateral triangle spans 
the two irreducible representations A and E, the latter 
one consisting of two associtcd one-dimensional represen- 
tations. To construct the SALC's, one simply applies the 
three projection operators V"^, , and to one of 
the original AO's, let's say the 0i, 

P^0i « (l)£0i + (l)(730i + (1)6-101 

= (1)01 + (1)02 + (1)03 

= 01 +02 + 03, (18) 

^^'01 ^ (l)£0i + (e)C'30i + (e*)C'|0i 

= 01 + £02 + £*03, (19) 

P^"01 « (l)£0i + (£*)(7301 + {s)Cl<Pi 

= 01+£*02+£03. (20) 

The A SALC in Eq. (18) has precisely the same form 
as the ^01 MO in Eq. (14), which was determined via a 
solution of the Hiickel equation (13). The two E SALC's 
[Eq. (19) and Eq. (20)], however, arc complex functions 
and do not coincide with the real ^02 and ^3 found above 
[Eq. (15) and Eq. (16)]. As we will see in section IV, these 
complex SALC's agree with BS UHF orbitals obtained in 
the case of an applied magnetic field. On the other hand, 
a set of two real and orthogonal SALC's that spans the 
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FIG. 5. The sS-UHF solution exhibiting breaking of the 
circular symmetry for A'' = 3 and Sz = 1/2 at Rw = 10 and 
-6 = 0. (a-b): real orbitals for the two spin-up electrons, (c): 
real orbital for the single spin-down electron, (d): total elec- 
tron density, (e): spin density (difference of the spin- up minus 
the spin-down partial electron densities). The choice of the 
remaining parameters is: htoo = 5 meV and m* = 0.067me. 
Distances are in nanometers. The real orbitals are in 10~^ 
nm"i and the densities (ED and SD) in 10"" nm"^. The 
arrows indicate the spin direction. 

E representation can be derived fron Eq. (19) and Eq. 
(20) by simply adding and substracting the two complex 
ones. This procedure recovers immediately the real V'2 
and discussed earlier. 

We stress here that the UHF orbitals of Fig. 4 are 
canonical (see section II. A) . As is well known from Quan- 
tum Chemistry,^^ in general the canonical spin orbitals 
will be spread out over the different sites (atoms) of a 
natural molecule and will form a basis for the irreducible 
representations of the symmetry group of the molecule. 
Once the canonical orbitals are available, there is an in- 
finite number of noncanonical spin orbitals that span 
reducible representions of the symmetry group of the 
molecule and can be obtained via a unitary transforma- 
tion of the canonical set. We remind the reader that 
noncanonical spin orbitals are solutions of a generalized 
HF equation involving off-diagonal elements Eij in the 
matrix formed out of the HF orbital energies (see section 



II. A). Naturally, the unitary transformation leaves the 
UHF determinant and total energy invariant. In partic- 
ular there is a unitary matrix that transforms the canon- 
ical spin orbitals to the fully localized AO's; i.e., for the 
iV = 3 and Sz = 3/2 case, such a unitary matrix trans- 
forms the canonical t/j^'s {i — 1, 2, 3) to the noncanonical 
AO's (j)iS. 

Note that such noncanonical orbitals for = 3 were 
recently used^^ to formulate a non-selfconsistent variant 
to the Pople-Nesbet HF equations listed in section II. A. 
This variant relied on the manifestation of spontaneous 
symmetry breaking that was discovered earlier via our 
selfconsistent UHF results. Notice, however, that Ref. 
51 obtained an incomplete wave function, since the com- 
panion step of restoration of the rotational symmetry was 
not considered (see section VI below). 

B. The = 1/2 partially spin polarized case 

In Fig. 5, we display the sS-UHF symmetry-violating 
orbitals (a— c) for the case of a partially polarized QD 
at Rw = 10 and B = Q with two spin-up and one spin- 
down electrons (A^ = 3 and = 1/2). The UHF orbital 
energies of these electrons are: (a) 45.350 and (b) 46.515 
meV for the spin-up orbitals, and (c) 45.926 meV for 
the spin-down orbital. An inspection of Fig. 5 reveals 
that these orbitals have retained the nodal structure of 
the corresponding independent-particle-model orbitals in 
the familiar Is^lpx configuration; namely, (a) and (c) are 
nodeless, while (b) exhibits a single nodal line. Apart 
of this property, however, and in consonance with the 
Sz = 3/2 case studied in a previous subsection, the BS 
UHF orbitals again differ drastically from the ones asso- 
ciated with the independent particle model; again they 
are associated with three sites within the QD arranged 
in an equilateral triangle. In contrast to the fully po- 
larized case, however, there are no linear combinations 
of atomic orbitals involving all three vertices of the tri- 
angle. Indeed the single spin-down electron remains by 
itself as an unmodified AO, while only the two spin-up 
electrons combine to form LCAO MO's. This behav- 
ior here is a special case of a general property of the 
sS-UHF, i.e., only AO's associated with the same spin 
direction can in principle combine to form LCAO MO's. 
This property, however, does not extend to the General- 
ized Hartree-Fock'^^'^^ which incorporates the additional 
unresctriction that the z-projection {Sz) of the total spin 
is not preserved and it is not a good quantum number. 
(Notice that unlike the practice in this paper, Ref. 31 
uses the term UHF for the Generalized HF.) 

In the same spirit with the treatment of the fully po- 
larized case in subsection HI. A, and taking into consider- 
ation the decoupling of the two different spin directions 
in the sS-UHF, one can write a corresponding Hiickel 
matrix equation for the = 3 and Sz — 1/2 case as 
follows. 
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FIG. 6. The S-UHF solution exhibiting breaking of the 
circular symmetry for A'' = 3 and = 3/2 at Rw = 10 
and B — 2 T. (a-c): orbitals (modulus square) for the three 
spin-up electrons, (d): total electron density. The choice of 
the remaining parameters is: huio = 5 meV, m* = 0.067me, 
and g* = —0.44. Distances are in nanometers and the densi- 
ties (orbital and ED) in 10~* nm~^. The arrows indicate the 
spin direction. 

( ~(3 f o] (h]=Elh]- (21) 

V \f,j \f,j 

From the eigenvalues and eigenvectors of (21), one 
finds the following three LCAO-MO's: (a) = (0i + 
<j>2)/V2 with energy Ei = e - f3, (b) ip2 = {(f>i - 4i2)/V^ 
with energy E2 = e + (3, and (c) -03 — 4>z with en- 
ergy i?3 — e. The structure of these three LCAO- 
MO's agrees very well with the corresponding symmetry- 
violating UHF orbitals for the = 3 and 6*2 = 1/2 case 
displayed in Figs. 5(a) — 5(c). This agreement extends 
also to the level diagrams of the corresponding orbital en- 
ergies. (Using the UHF values for Ei and E2, one finds 
(3 = 0.582 meV and e = 45.932 meV « £^3.) 

Concerning the underlying symmetry-group structure 
of the UHF orbitals in Fig. 5, we observe that (unlike the 
fully polarized case) the two rotations, i.e., C3 (rotation 
by 27r/3) and C| (rotation by 47r/3), are not part of the 
symmetry operations of the relevant point group [due to 
the dissimilarity between spin- up and spin-down orbitals; 
this can be seen clearly through inspection of the spin 
density in Fig. 5(c)]. Because of the 2D character of the 
dot, the only symmetry operations for the A^ = 3 and 
Sz = 1/2 case are the identity E and a single reflection, 
al, through the vertical plane passing through the spin- 
down electron. Such a group {E, a^} is a subgroup of the 
familiar C2V group associated with the 3D H2O molecule 
(observe that the O atom corresponds to the spin-down 
electron and the two H atoms correspond to the two spin- 
up electrons). According to Ref. 33 (see p. 181), the 



representation F formed by the three original AO's can 
be reduced to irreducible A2 and Bi representations as 
T ^ A2 + 2Bi. By applying projection operators [see 
Eq.(17)] to the 01 AO and using the character Table H, 
one finds the following two normalized SALC's, 

V'A. = (01 - 02)/\/2, (22) 

and 

i^B, = (01 + 02)/\/2. (23) 

Naturally the second SALC of Bi symmetry is 

V'k = 03. (24) 

Once more, we stress the fact that the SALC's [Eqs. 
(22)— (24)] derived above via symmetry-group theory 
have the same structure as the canonical UHF orbitals 
displayed in Fig. 5. Observe further that in the General- 
ized HF (as also in the case of the H2O molecule and the 
allyl anion) the two SALC's of Bi symmetry are allowed 
to couple, producing more complicated orbitals and en- 
ergy level diagrams (see p. 1037 of Rcf. 31 and p. 183 in 
Ref. 33) 



IV. THREE ELECTRONS IN A FINITE 
MAGNETIC FIELD 

A. The Sz =3/2 fully polarized case 

In this section, we study the case of three fully po- 
larized electrons under a magnetic field B. Since for 
B ^ 0, the BS UHF orbitals are necessarily complex 
functions. Fig. 6 displays the modulus square of these 
orbitals. The UHF total ED displayed in Fig. 6(d) 
and the modulus square of the orbitals exhibit an ap- 
parent Csy symmetry as was the case at -B = (sec- 
tion HLA). However, the phases of the complex orbitals 
at B ^ contribute to a modification of the symme- 
try group. This modification has been studied earlier 
for the case of infinite crystalline systems with periodic 
space lattices where the electrons occupy Bloch orbitals, 
and such studies have led to the consideration of two 
physically equivalent group structures, namely the ray 
groups^^ and the magnetic translation groups. In our 
case of a finite periodic crystallite, corresponding mag- 
netic rotation groups would be straightforward to con- 
sider. However, in order to appreciate the modifications 
introduced by the magnetic field, it will be simpler to 



TABLE II. Character table for the group {E, ai} 
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FIG. 7. The S-UHF orbital energies (in meV) for 
N = S,Sz = 3/2 and Rw = 10 as a function of the magnetic 
field B (in Tesla) , exhibiting a prominent Aharonov-Bohm os- 
cillation. The choice of the remaining parameters is: ?uJo = 5 
meV, m* = 0.067me, and g* = -0.44. 



modify the Hiickel (tight-binding) hamiltonian according 

to the Harpcr-Pcicrls prescription,"'''''''''*' which accounts 
for the magnetic gauge transformation when moving from 
one crystalHne site to another. 

Thus according to Pcierls and Harper, the proper 
atomic orbitals (j)j's for B ^ (centered at Rj) are the 
real (f)/s multipUed by an appropriate phase as follows, 



,(r;R,) = <^,(r;R,)et|A(R. 



(25) 



Because of this position-dependent phase in the AO's, the 
hopping matrix elements Hij (see section III.A) are now 
complex, and the Hiickel equation (13) for three electrons 
is modified as follows: 



e 



-/3e'' 
-/Je- 




(26) 



where O = flij = (e/fi.c)A(R.j — Ri)'Ri, = 
(1,2), (2,3), (3,1), with R^, fc = 1,2,3 being the posi- 
tions of the vertices of the equilateral triangle. Notice 
that n = (27r$)/(3$o); where $ is the total magnetic 
flux through the equilateral triangle and $o = hc/e is 
the unit flux. 

From the eigenvectors of equation (26), one finds the 
following LCAO-MO's, 



■t/jl OC + (j)2 + 03, 



>3, 



(27) 



(28) 



(29) 



with corresponding orbital energies. 
El = e - 2/JcosO, 



2/3cos[(27r/3-hO)], 



(30) 
(31) 



iJs =e-2/3cos[(27r/3-ri)]. (32) 

Substituting the specific value for Cl given above, one 
can write the eigenvalues (30-32) in a more symmetric 
compact form. 



Ei 



2/3 cos 



27r,. 



,i = 1,2,3. (33) 



Since the original AO's do not practically overlap for 
Rw = 10, the phases in front of the (j)i's in Eqs. (27-29) 
do not contribute substantially to the modulus square of 
the orbitals. As a result, for all values of B, all three 
orbitals exhibit similar orbital densities that are approx- 
imately equal to + 02 + (pl- Observe that this agrees 
very well with the behavior of the canonical UHF orbitals 
(modulus square) at i? = 2 T displayed in Fig. 6. At zero 
magnetic field B = 0, the LCAO-MO's in equations (27- 
29) reduce to the specific form given earlier in Eqs. (18- 
20) of section III.A. We stress here that in section III.A 
these LCAO-MO's were derived from arguments based 
exclusively on the group theoretical structure of the C^y 
symmetry group. 

Naturally, when = 0, the orbital energies in Eq. (33) 
reduce to the corresponding B = result derived in sec- 
tion III.A, namely Ei = e - 20 and E2 = E3 = e + (3. 
Notice, however, that for arbitrary values of B, the de- 
generacy between E2 and E3 is lifted. In addition, the 
three energies Ei, E2, and E3 in Eq. (33) exhibit promi- 
nent Aharonov-Bohm (AB) oscillations. It is interesting 
to compare this behavior to the behavior of the calculated 
canonical UHF orbital energies for Rw = 10. These UHF 
orbital energies as a function of B are displayed in Fig. 
7. An inspection of Fig. 7 reveals that the UHF orbital 
energies do exhibit (as expected) an Aharonov-Bohm os- 
cillation as a function of B. However, these oscillations 
are more complicated from what can be simply antici- 
pated from the analytic formulas in Eq. (33). Namely, 
the amplitude of the UHF AB oscillations decreases as B 
increases. This behavior is due to a decrease in the hop- 
ping parameter /3, which results from the spatial shrink- 
age of the Gaussian-type UHF orbitals as a function of 
B. Eventually, for B — > 00, all three energies are degen- 
erate. We notice that complete degeneracy of all UHF 
orbitals for any A'' appears also in the B = 0, Rw — * 00 
limit. 

The electronic structure of the UHF fully-polarized 
three-electron molecule in a magnetic field, which was 
discussed above and which exhibits Aharonov-Bohm os- 
cillations, does not have an analog in the realm of natural 
molecules. However, apart from the B-dependence of /?, 
it agrees in a remarkable way with the "noninteracting 
spectra" of the artificial molecules that can be formed 
out of ID ring arrays of single QD's.^^ 
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FIG. 8. The sS-UHF solution exhibiting breaking of the cir- 
cular symmetry for iV = 3,3;, = 1/2 at Rw = 10 and B = 2 T. 
(a-b): orbitals (modulus square) for the two spin-up electrons, 
(c): orbital (modulus square) for the spin-down electron, (d): 
total electron density. The choice of the remaining parameters 
is: huJo = 5 meV, m* — 0.067me, and g* = —0.44. Distances 
are in nanometers and the densities (orbital and ED) in 10"* 
nm~'^. The arrows indicate the spin direction. 



B. The Sz = 1/2 partially polarized case 

Fig. 8 displays the BS UHF orbitals (modulus square) 
and the total ED for the partially polarized N — 3, — 
1/2 case in a magnetic field i? = 2 T. As with the B = Q 
case (section III.B), the spin-down orbital is decoupled 
from the two spin-up ones. As a result the corresponding 
Hiickel matrix equation is of the form 



e 

-pe- 




in 














(34) 









where in general e' 7^ e due to the energy difference be- 
tween the two spin directions introduced by the Zeeman 
term. From the solutions of Eq. (34), one finds the fol- 
lowing LCAO-MO's: (a) ^1 = (e*"0ij- h)/V2 with 
energy Ei = e — /3, (b) V'2 = (e'^(/)i — (/)2)/V2 with en- 
ergy E2 = e + f3, and (c) ^3 — 03 with E3 = e'. The 
total electron density constructed out of these LCAO- 
MO's is again of the form (f)1 + (j)2 + 03 (compare with 
section IV. A). The corresponding UHF orbitals (mod- 
ulus square) and ED displayed in Fig. 8 are obviously 
conforming to these forms. 

Concerning the Hiickel orbital energies Ej, j — 1, 2, 3, 
we note that they do not depend on the magnetic field 
B through ft. As a result, unlike the previous case of 



the fully polarized electrons, AB oscillations should not 
develop in the UHF orbital energies. To check this pre- 
diction, we dispaly in Fig. 9 the UHF orbital energies 
as a function of B. In contrast to Fig. 7, AB oscilla- 
tions are absent in Fig. 9, a behavior which apparently 
relates to the fact that no UHF orbital covers the area 
of the equilateral triangle (the single spin-down orbital 
does not couple to the two spin- up ones, which lie on a 
straight hue). 



V. SIX ELECTRONS AT ZERO MAGNETIC 
FIELD 

We discuss now the case of six fully polarized electrons 
in zero magnetic field. The corresponding total S-UHF 
electron density for Rw = 15 (k = 1.2730) is displayed 
in Fig. 10(a) (bottom frame). Unlike the case of smaller 
numbers of particles with iV < 5, six electrons is the 
smallest system that forms a Wigner molecule with a two- 
ring arrangement. Such a ring arrangement is denoted 
by (1, 5) to distinguish it from a single-ring arrangement 
(0,iV). 

Naturally, single-ring molecular arrangements (0, N) 
are familiar from the Quantum Chemistry of carbocyclic 
systems (Ref. 33, see also sections III and IV). The more 
complicated (1,5) arrangement, however, is a molecular 
structure unknown to traditional Chemistry. Neverthe- 
less, contact to Organic Chemistry can be retained by 
observing that the (0,5) outer ring has a group symme- 
try similar to the cyclic hydrocarbon C^H^. As a result, 
and in direct analogy with the molecule (see p. 152 

in Ref. 33), the SALC's of the (0,5) arrangement are as 
follows, 



+ 02 + 03 + 04 + 05), 



(35) 
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FIG. 10. The canonical S-UHF real orbitals for TV = 6 and 
~ 3, and for Rw = 15 and B = 0. (a): the total electron 
density, (b-c), middle row: the two orbitals of A symmetry, 
(d-e), top row: the two degenerate orbitals of E2 symmetry. 
The choice of the remaining parameters is: huo — 5 meV 
and m* — 0.067me. Distances are in nanometers. The real 
orbitals are in 10~^ nm~^ and the electron density in 10"* 
nm^^. The arrows indicate the spin direction. 



r (j>2 COS 6 - 

. COS 26 + t 



1 COS 26 



05 cos f 




ti2 sin 9 + (j>3 sin 26 
- (j)4 sin 26 — 05 sin 0) , 



1+02 cos 26 + 03 cos ( 



cos6' + 05 cos2(?), 



ip{E2b) = (02 sin 26* - 03 sin 6* 

+04 sin6' — 05 sin 20), 



(36) 



(37) 



(38) 



(39) 



where 6 = 2tt/5. The corresponding orbital energies are 
e — 2/3 for the single orbital of A symmetry, e — (2 cos 9)/3 
for the two degenerate orbitals of Ei symmetry, and e — 
(2cos26')/3 for the remaining two degenerate orbitals of 
E2 symmetry. 

Returning back to the case of the (1,5) ring arrange- 
ment, we notice that the sixth AO, 06, at the center is of 
A symmetry, and thus it can only couple to a MO of the 
(0,5) ring with the same symmetry, namely the orbital 
-ipiA) in Eq. (35). As a result, both the (0,5) and the 
(1,5) ring arrangements share the same four MO's of Ei 
and E2 symmetry. 

The coupling matrix element between the 06 and 'ip{A) 
orbitals is given by, 

5 

J 06i?^(A)dr = -L ^ y 06H0fcrfr 



-V5d. 



(40) 



To find the MO's of the (1,5) ring with A symmetry, 
we need to solve the 2x2 matrix equation. 



e-2/3 -V56 
e 



91 

92 



E 



91 

92 



(41) 



We note that, due to the different coordination and 
distances, the quantities 6 and associated with the cen- 
tral AO are different from the corresponding quantities 
/3 and e associated with the AO's of the outer ring. 

Using the notation. 



v/20(52 + (?- 6 + 2/3)2, 



(42) 



the two solutions of the matrix equation (41) have ener- 
gies (?+ e — 2/3 + Q)/2 and eigenvectors (unnormalized) 
{e — e + 2f3 ± Q, ^/205}, respectively. Accordingly, one 
MO of the (1,5) ring is constructed by adding and the 
other by subtracting a fraction of the ip{A) orbital from 
the central AO. This behavior agrees very well with the 
two S-UHF orbitals displayed in the second row of Fig. 
10 [(b) and (c)]. 

Concerning the S-UHF orbitals displayed in the third 
row of Fig. 10 [(d) and (e)], we notice that they are de- 
generate in energy and that they agree very well with the 
two MO's of E2 symmetry. Indeed the S-UHF orbital in 
Fig. 10(d) exhibits five total humps, three of them posi- 
tive and the other two negative, in remarkable agreement 
(apart from an overall sign) with the MO in Eq. (38) [note 
that cos 6* = 0.3090 > and cos 26* = -0.8090 < 0]. In 
particular, counterclockwise, the polarities of the humps 
in Fig. 10(d) are (—,+,—,—,+), differing only by an 
overall sign from the corresponding polarities of the MO 
in Eq. (38). Even more impressive is the fact that there 
is quantitative agreement regarding the absolute heights 
of the humps in these two orbitals (see the values of 
cos 6 and cos26' listed above). The other S-UHF or- 
bital in Fig. 10(e), exhibits a total of four humps, two 
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of them positive and the other two negative, and hav- 
ing an alternating (+,—,+,—) arrangement. This is 
again in remarkable agreement with the second MO of 
E2 symmetry in Eq. (39), since sin0 > and sin 20 > 0. 
Additionally, we note that the agreement between the 
UHF orbital in Fig. 10(e) and the MO in Eq. (39) ex- 
tends further to the absolute heights of the humps, since 
sin 6* = 0.9511 > 0.5878 = sin 261. 

Finally, there arc two other degenerate UHF orbitals 
that are not displayed in Fig. 10. They are not identical 
to the V(-Eia) and ip{Eib) SALC's in Eqs. (36) and (37), 
but we have checked that they span the Ei irreducible 
representation. 

VI. RESTORATION OF CIRCULAR SYMMETRY 
AND EXACT SPECTRA 

A. Group structure and sequences of magic angular 
momenta 

In the previous sections, we demonstrated that the BS 
UHF determinants and orbitals describe indeed 2D elec- 
tronic molecular stuctures (Wigner molecules) in close 
analogy with the case of natural 3D molecules. How- 
ever, the study of the WM's at the UHF level restricts 
their description to the intrinsic (nonrotating) frame of 
reference. Motivated by the case of natural atoms, one 
can take a subsequent step and address the properties 
of collectively rotating WM's in the laboratory frame 
of reference. As is well known, for natural atoms, this 
step is achieved by writing the total wave function of the 
molecule as the product of the electronic and ionic par- 
tial wave functions. In the case of the purely electronic 
WM's, however, such a product wave function requires 
the assumption of complete decoupling between intrin- 
sic and collective degrees of freedom, an assumption that 
might be justifiable in limiting cases only. 

As we demonstrated earlier, in the framework 
of the BS UHF solutions this companion step can be 
performed by using the post-Hartrec-Fock method of 
restoration of broken symmetries'^^ (RBS) via projec- 
tion techniques (PT's). Examples demonstrating the 
RBS method have been presented by us in two cases: 
(I) The ground state (with angular momentum 7 = 0) 
of two interacting electrons in a parabolic quantum dot 
in the absence of a magnetic field; (II) The yrast ro- 
tational band (sec section I.C and precise definition in 
Ref. 58) of a system of N interacting electrons in high 
magnetic fields^^'^^ (fractional-quantum-Hall regime). In 
both cases, we showed that the RBS method (as adapted 
to the case of 2D BS UHF solutions) yields correlated 
(multideterminantal) many-body wave functions that ap- 
proximate very well the corresponding exact solutions. 
In particular, in the latter case, our use^^'^"^ of the RBS 
method yielded analytic expressions for the correlated 
wave functions that offer a better description of the A''- 



electron problem in high B compared to the Jastrow- 
Laughlin'"' expression. 

In this section, we will not proceed any further with ex- 
plicit numerical or analytic derivations of additional RBS 
wave functions. Instead, we will use the RBS approach to 
illustrate through a couple of concrete examples how cer- 
tain universal properties of the exact solutions, i.e., the 
appearance of magic angular momenta in the exact rota- 
tional spectra, ^" relate to the symmetry broken UHF 
solutions. Indeed, we will demonstrate that the magic an- 
gular momenta are a direct consequence of the symmetry 
breaking at the UHF level a,nd that they are determined 
fully by the molecular symmetries of the UHF determi- 
nant.^^ 

As an illustrative example, we have chosen the rela- 
tively simple, but non trivial case, of A'' = 3 electrons. 
For 5 = 0, both the 5*^ 1/2 and 5^ = 3/2 polariza- 
tions can be considered. We start with the Sz = 1/2 
polarization, whose BS UHF solution (let's denote it by 
I []])) was presented in section III.B and which exhibits 
a breaking of the total spin symmetry in addition to the 
rotational symmetry. We first proceed with the restora- 
tion of the total spin by noticing that | []]) has a point- 
group symmetry lower (see section III.B) than the C^v 
symmetry of an equilateral triangle. The symmetry, 
however, can be readily restored by applying the projec- 
tion operator (17) to | J,tT) and by using the character 
table of the cyclic C3 group (see Table I). Then for the 
intrinsic part of the many-body wave function, one finds 
two different three-determinantal combinations, namely 

<M = I itt) + e'^^^/'l Tit) + e-^-^'^\ TTi), (43) 
and 

*Str(7o) = I itt) + e-2"'/'| Tit) + e'"'/'| TTi), (44) 

where 70 = denotes the azimuthal angle of the vertex 
associated with the original spin-down orbital in | J,tT)- 
We note that the intrinsic wave functions ^j^tj. and <I>i^tr 
are eigenstates of the square of the total spin operator 
(S = X]?=i ^i) with quantum number s = 1/2. This 
can be verified directly by applying to them.^° 

To restore the circular symmetry in the case of 
a (0,N) ring arrangement, one applies the projection 
operator, ^^'^^ 

/■27r 

2'kVi= / rf7exp[-i7(L-7)] , (45) 
Jo 

where L = X^^Li h is the operator for the total angular 
momentum. Notice that the operator Vi is a direct gen- 
eralization of the projection operator (17) to the case of 
the continuous cyclic group C^o [the phases exp(i7/) are 
the characters of Coo]- 

The RBS projected wave fimction, 'I'rbSj (having both 
good total spin and angular momentum quantum num- 
bers) is of the form, 
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27r*RBS = / rf7^.^tr(7)e'^^ (46) 
Jo 

where now the intrinsic wave function [given by Eq. (43) 
or Eq. 44)] has an arbitrary azimuthal orientation 7. We 

note that, unhkc the phcnomcnological Eckardt- frame 
model^^'^^ where only a single product term is involved, 
the RBS wave function in Eq. (46) is an average over all 
azimuthal directions of an infinite set of product terms. 
These terms are formed by multiplying the UHF intrin- 
sic part $^^^^(7) with the external rotational wave func- 
tion exp(i7/) (the latter is properly characterized as "ex- 
ternal" , since it is an eigenfunction of the total angular 
momentum L and depends exclusively on the azimuthal 
coordinate 7). 

The operator _R(27r/3) = cxp(— j27rL/3) can be applied 
onto ^'rbs ill two different ways, namely cither on the in- 
trinsic part or the external part cxp(z7/). Using Eq. 
(43) and the property .R(27r/3)<I>j^i^ = exp(-27ri/3)$j^t^, 
one finds, 

^(27r/3)*RBs = exp(-27ri/3)*RBS, (47) 

from the first alternative, and 

i?(27r/3)*RBs = exp(-27r/V3)*RBS, (48) 

from the second alternative. Now if 5'rbs 7^ 0, the only 
way that Eqs. (47) and (48) can be simultaneously true 
is if the condition exp[27r(/ — 1)V3] = 1 is fulfilled. This 
leads to a first sequence of magic angular momenta asso- 
ciated with total spin s = 1/2, i.e., 

/ = 3fc + 1, fc = 0, ±1, ±2, ±3, ... (49) 

Using Eq. (44) for the intrinsic wave function, and fol- 
lowing similar steps, one can derive a second sequence of 
magic angular momenta associated with good total spin 
8 = 1/2, i.e., 

7 = 3A; - 1, A; = 0, ±1, ±2, ±3, ... (50) 

In the fully polarized case, the UHF determinant was 
described in section III. A. This UHF determinant, which 
we denote as | ttt )> is already an eigenstate of 
with quantum number s — 3/2. Thus only the rota- 
tional symmetry needs to be restored, that is, the intrin- 
sic wave function is simply $(^^^(70) = | TTT )• Since 
i?(27r/3)<I'j^jj. = $|^iti.j the condition for the allowed an- 
gular momenta is exp[— 27r/i/3] = 1, which yields the 
following magic angular momenta, 

7 = 3fc, fc = 0,±l,±2,±3, ... (51) 

We note that in high magnetic fields only the fully po- 
larized case is relevant and that only angular momenta 
with fc > enter in Eq. (51) (see Ref. 16). In this case, 
in the thermodynamic limit, the partial sequence with 
fc = 2^ -h 1, g' = 0, 1, 2, 3, ... is directly related to the odd 



TABLE III. Case of A'^ = 6 electrons in high magnetic field 
B: Total interaction energies in the lowest Landau level of 
REM and oxact-diagonalization wave functions for various 
magic angular momenta I of the yrast band. The REM func- 
tions are analytically specified RBS wave functions derived in 
Ref. 16. The percentages within parenthesis indicate relative 
errors. Energies in units of /kIb, where k is the dielec- 
tric constant and Ib = yJhc/eB is the magnetic length. For 
details concerning the exact-diagonalization method and the 
REM wave functions, see Ref. 17. For additional values of 7, 
see Ref. 17. 

7 REM EXACT 

70 2.3019 (0.85%) 2.2824 

80 2.1455 (0.71%) 2.1304 

90 2.0174 (0.60%) 2.0053 

100 1.9098 (0.51%) 1.9001 

110 1.8179 (0.45%) 1.8098 

120 1.7382 (0.40%) 1.7312 

130 1.6681 (0.36%) 1.6621 



filling factors u = l/(2g -|- 1) of the fractional quantum 
Hall effect [via the relation v = N{N - l)/(27)]. This 
suggests that the observed hierarchy of fractional filling 
factors in the quantum Hall effect may be viewed as a 
signature originating from the point group symmetries 
of the intrinsic wave function $intr, and thus it is a man- 
ifestation of symmetry breaking at the UHF mean-field 
level. 



B. Quantitative description of the yrast band 

The usefulness of the RBS wave functions [Eq. (46)] 
is not limited to deriving universal properties of the ex- 
act spectra, like the sequences of magic angular momenta 
[sec section VI. A]. As we demonstrated in earlier publi- 
cations, in the the regime of strong correlations, the RBS 
wave functions approximate very well the corresponding 
exact many-body wave functions. 

Indeed, in Ref. 15 we offered (as a function of Rw) 
a systematic comparison between the RBS and exact 
ground-state (7 = 0) energies at i? = for A'' = 2 elec- 
trons in a parabolic QD. For Rw = 19.09, we found that 
the relative error was approximately 0.7%. Furthermore 
in Ref. 16, for the case of high B, we derived analytic 
RBS wave functions, named "REM wave functions". As 
we showcd"'^^ explicitly for the case of A^ = 6 electrons, 
the radial electron densities associated with the REM 
functions accurately reproduce the ones extracted from 
exact-diagonalization calculations. 

In this subsection, we offer additional examples per- 
taining to the ability of the RBS wave functions to repro- 
duce the exact yrast spectra of parabolic QD's. In par- 
ticular, Table III lists the REM and exact yrast energies 
in the range of magic angular momenta 70 < 7 < 130 for 
N = 6 electrons in high B. Details concerning the REM 
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TABLE IV. Case of iV = 2 electrons in a parabolic QD 
at B = : Total energies of RBS and exact wave functions 

for various magic angular momenta / of the yrast band. The 
percentages within parenthesis indicate relative errors. The 
choice of remaining parameters is: fujo = 5 meV, k = 1 
{Rw = 19.09), m* = 0.067me. Energies in units of meV. For 
details concerning the method for finding exact solutions to 
the two-electron problem, see Ref. 14. For details concerning 
the calculation of the RBS yrast spectrum, see Ref. 15 and 
the Appendix. 



I RBS EXACT 

52.224 (0.75%) 51.831 

1 52.696 (0.77%) 52.292 

2 54.086 (0.88%) 53.615 

3 56.240 (1.04%) 55.654 

4 59.065 (1.39%) 58.255 

5 62.065 (1.27%) 61.285 

6 63.911 (1.13%) 64.642 



wave functions and the exact-diagonalization method in 
the lowest Landau level are given in Ref. 17, and they 
will not be repeated here. 

The RBS and exact yrast spectra (0 < / < 6) for the 
case of A'' = 2 electrons at B = and Rw = 19.09 are 
given in Table IV. Details concerning their calculation are 
given in Ref. 15 and in the Appendix (where we present 
the final formula for calculating RBS energies for both 
even and odd angular momenta). 

We note that the relative errors in both Table III and 
Table IV are small (smaller than 1% in the majority of 
cases). 

VII. SUMMARY 

In this paper, we have introduced a group theoreti- 
cal analysis of broken-symmetry UHF orbitals and total 
electron densities in the case of single 2D semiconductor 
QD's. This analysis provided further support for our ear- 
lier interpretation^'"'^'''"'^^'"'^^ concerning the spontaneous 
formation of collectively rotating electron (or Wigner) 
molecules. Indeed the group-theoretical analysis enabled 
us to unveil further deeper analogies between the elec- 
tronic structure of the Wignc;r molecules and that of the 
natural 3D molecules. In particular these deeper analo- 
gies are: (I) The breaking of rotational symmetry results 
in canonical UHF orbitals that are associated with the 
eigenvectors of a molecular-type Hiickel hamiltonian with 
sites at positions specified by the equilibrium configura- 
tion of the classical A''-electron problem; (II) The broken- 
symmetry canonical UHF orbitals transform according to 
the irreducible representations of the point group speci- 
fied by the discrete symmetries of this classical molecular 
configuration; (III) The WM's formed out of the broken- 
symmetry UHF solutions can rotate, and the restoration 
of the total-spin and rotational symmetries results (in ad- 



dition to the ground state) in states defining the lowest 
rotational bands (i.e., yrast bands) of the WM's; (IV) 
The breaking of the circular symmetry results in lower- 
ing of the symmetry. This is expressed by the discrete 
point-group symmetry of the UHF wave function and it 
underlies the appearance of sequences of magic angular 
momenta (familiar from exact-diagonalization studies) in 
the excitation spectra of single QD's. 

Since exact-diagonalization methods arc typically re- 
stricted to small sizes with N < 10, the two-step method 
of breakage and subsequent restoration of symmetries 
offers a promising new venue for accurately describing 
larger 2D electronic systems. A concrete example of 
the potential of this approach is provided by Ref. 16, 
where our use of the the symmetry-breaking/symmetry- 
restoration method yielded analytic expressions for cor- 
related wave functions that offer a better description of 
the A^-electron problem in high magnetic fields compared 
to the Jastrow-Lauglilin'' expression. 

Furthermore, the group-theoretical analysis strongly 
suggests an interesting simplified variant approach for 
carrying out the first step of symmetry breaking. This 
variant rests on the observation that, in all cases of 
WM's, the broken-symmetry UHF orbitals arc generic 
linear combinations of Gaussian-type functions [with a 
proper phase for B 0, see Eq. (25)] specified simply 
by their width a and positions R^'s from the center of 
the QD. The linear combinations can be fuUy specified 
from the group theoretical analysis of the appropriate 
classical equilibrium configurations,^^ and a determinant 
of the corresponding LC AO-MO 's can readily be writ- 
ten down. Then a simple variational calculation of the 
minimum total energy of this determinant will yield the 
parameters a and Rj 's without the need to carry out the 
self-consistent UHF iterations. This simplified approach 
could treat even larger sizes without major loss of accu- 
racy. Added accuracy can then be obtained through the 
subsequent step of restoration of the broken symmetries. 
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APPENDIX A 

For the case of iV = 2 electrons in a parabolic QD at 
B = 0, we reported in Ref. 15 the RBS formulas for cal- 
culating energies of yrast-band states with even angular 
momenta /. These formulas [see Eqs. (11)-(13) in Ref. 
15] were generated via a projection of the "singlet" UHF 
determinant. The corresponding RBS formulas for odd 
values of / are generated via a projection of the triplet 
UHF state. 
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In this Appendix, wc present the formulas covering 
both even and odd angular momenta. They are: 

/■27r / f2TT 

Srbs(/) = h{^)e'^U^ n(7)e*T^rf7, (Al) 



with 



(A2) 



and 



n(7) = SusSvt ± SutSy 



(A3) 



where the upper signs apply in the case of even /'s and 
the lower signs in the case of odd I's. s(r) and t{r) are 
the initial u{r) and v{r) broken-symmetry UHF orbitals 
rotated by an angle 7, respectively. Vuhsl and Vuyts sue 
two-body matrix elements of the Coulomb repulsion, and 
Sus, etc., are the overlap intergrals. 
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